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ABSTRACT. Working directly from the 3D magnetohydrodynamical equations and entirely in physi- 
cal scales a scenario is identified wherein the kinetic and magnetic enstrophies in a turbulent plasma 
exhibit a cascade over a range of scales the lower limit of which is a modified Kraichnan-type scale. 
Sufficient conditions include a hybrid geometric/ smoothness assumption on the vorticity, a restraint 
on the scales over which the cascade transpires, and technical localization and modulation require- 
ments on the evolution of the vorticity and current. Scale locality of the enstrophy flux is also exam- 
ined. 

. i. Introduction. 

Q-r 

In a series of recent publications, R. Dascaliuc and Z. Grujic established the existence of various 
turbulent cascades as well as the locality of the corresponding fluxes working in physical scales of 
incompressible fluid flows. The results rigorously follow directly from the Navier-Stokes equa- 
tions (cf. MDaGrlL MDaGr2L ||DaGr3L and !IDaGr4ll ). This was achieved via a novel process of 
ensemble averaging which acts as a detector of significant sign-fluctuations of a physical density of 
interest. In the context of turbulent cascades, the a priori sign-changing density of interest is flux; 
if a quantity such as energy or enstrophy is being transferred from larger to smaller scales then 
the (suitably averaged) inward directed flux through shells of a given radius should be positive. 
That the setting is entirely physical provides a natural way to incorporate geometric information 
about the flow into the analytic framework. This was particularly useful in HDaGr3| . 

Because of the formal similarities between the two systems, mathematical results in the theory of 
NSE often illuminate and direct the study of MHD. For general remarks concerning the mathemat- 
ical theory of MHD see IISeTel . In this and a complementary paper, IBrGrlL the authors extend 
the framework of Dascaliuc and Grujic to study various cascades present in magnetohydrody- 
namic turbulence (see ]DB| for the physical theory of MHD turbulence). The focus of IBrGrl| are 
turbulent events at the u-b level. There, the total energy flux is found to cascade over a range of 
scales the lower bound of which is a modified Taylor microscale. Further results are given for 
energies transported by specific mechanisms (advection by the magnetic field, for instance). Re- 
garding the cascade for the fluid-advected kinetic energy plus fluid pressure flux, a role is identi- 
fied for a numerically and phenomenologically motivated scale dependent dynamic alignment (cf . 
IB051 IPeBolOl IPoBhlOl IGS951 ) which serves to partially deplete the non-linear coupling between 
the magnetic and velocity fields. In addition to studying turbulent cascades, HBrGrlH includes a 
scenario in which the predominant energy transfer is from the velocity to the magnetic field, re- 
sults regarding locality of the fluxes, and a note detailing implications for non-decaying turbulence 
in flows governed by 3D NSE. 

The present work identifies a set of physically reasonable assumptions which guarantee the exis- 
tence of magnetic and kinetic enstrophy cascades across a range of physical scales as well as the 
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scale locality between scales of these cascades. The framework and trajectory of our efforts fol- 
low the work presented in BDaGr3l but, due to the numerous non-linear coupled terms evident in 
the current-vorticity formulation of 3D MHD, significant technical adaptations are needed. As in 
HDaGr3l , a significant mathematical hurdle is to determine a range of physical scales over which 
the contribution of non flux-type terms to the ensemble averages can be interpolated between cer- 
tain integral-scale quantities in order to ensure their detractions are small. In the incompressible 
fluid case, the only term requiring such treatment is the vortex stretching term and it is here that 
the geometry of coherent vortex structures plays an important role. The 3D MHD setting is con- 
siderably more complicated. The current-vorticity formulation includes a number of critical order 
advective terms which are structurally similar to the vortex stretching term but also includes a 
constructurally unique term, Vui x Vbi . Treatment of this last term necessitates several modifica- 
tions to the assumptions prescribed in IIDaGr3L 

The most significant change between the assumptions in MDaGr3ll and our efforts here is necessi- 
tated because the purely geometric criteria motivated by the existence of coherent vortex struc- 
tures is not apparently helpful in the treatment of Vii; x Vfy. Modified assumptions have been 
detailed in MWuOOl for ideal MHD and in MHeXinll in order to establish comparable regularity re- 
sults as MCoFe93M and MBeBeM where the methods applied in MDaGr3l were pioneered. Our assump- 
tion is the hybrid geometric/ smoothness condition of MHeXinll which, instead of asserting local 
coherence on the directions of vorticity, requires that, for a point x lying in a region of high kinetic 
turbulence - i.e. | Vu| large - and for \y\ less than some critical scale, we have the inequality: 

\oj(x + y) — w{x)\ < \oj(x + y)\\y\ 5 . 

The above relationship, due to the fact that the strain tensor can be decomposed into singular 
integrals with kernels exhibiting certain cancellation properties over the unit sphere (cf. IICo94ID , 
results in the same localized depletion of nonlinearities present in MDaGr3ll for the vortex stretching 
term. It is worth highlighting that the hybrid geometric / smoothness condition is imposed solely 
on the vorticity field. 

The differences between the fluid and plasma cases are also evidenced by altered dependencies of 
the Kraichnan-type scale. From our treatment of Vui x Vfy we accumulate a dependence on an 
integral scale and energy-level space-time averaged quantity in addition to one at the enstrophy- 
level (the latter of these is the only one present in the fluid case). This is necessary to accommodate 
the local discrepancy between j and Vb. That is, although ||j||2 = ||V6||2, the same cannot be said 
for | \cj) j\ 1 2 and | \4> Vb\ I2 where is a smooth cut-off function having bounded support. 

This paper is organized as follows. In Section 2 we recall the theoretical apparatus by which we 
study turbulence cascades in physical scales. In Section 3 we contextualize the general discussion 
of Section 2 to the case of 3D MHD and present the main theorem, Theorem [5] Section 4 is con- 
cerned with establishing necessary bounds for the proof of Theorem|5] while Section 5 contains the 
actual proof. Finally, in Section 6, we include corollary results establishing the scale locality of the 
enstrophy fluxes. 

2. (K 1 ,K 2 ) -Covers and Ensemble Averages. 

The main purpose of this section is to describe how ensemble averaging with respect to (K\, K2)- 
covers of an integral domain B(0, Rq) can be used to establish essential positivity of an a priori sign- 
varying density over a range of physical scales associated with the integral domain. The appli- 
cation to turbulence is establishing the positivity of certain inward directed flux densities - i.e. 
the cascade is uni-directional from larger to smaller scales - as well as the near-constancy of the 
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averaged densities - i.e. the space-time averages over cover elements are all mutually comparable 
- across a range of scales. 

The ensemble averages will be taken over collections of spatio-temporal averages of physical den- 
sities localized to cover elements of a particular type of covering - a so called (K\, K-i)-cover - 
where the cover is over the region of turbulent activity. For simplicity this region will be taken as 
a ball of radius Rq centered at the origin and is henceforth referred to as the integral domain. The 
(Ki, i^-covers are now defined. 

Definition 1. Let K\, K2 £ N and < R < Rq. The cover of the integral domain B(0, Rq) by the n balls, 
{B(xi, R)}f = i is a (Ki, i^ 2 )-cover at scale R if, 

and, for any x £ B(0, Rq), x is contained in at most K~2 balls from the cover. 



In the hereafter all covers are understood to be {K\,K~2)- covers at scale R. The positive integers 
K\ and K~2 represent the maximum allowed global and local multiplicities, respectively. 

In order to localize a physical density to a cover element we incorporate certain refined cut-off 
functions. For the cover element around the point Xj, let 4>i{x, t) = rj(t)tp(x) where 77 £ C°°(0, T) 
and i\) £ Cg°(B(xi,2R)) satisfy 

(1) < 77 < 1, 77 = Oon (0,T/3), 7] = 1 on (2T/3,T), Mr" < 



r] 5 ~ T 



and, 



(2) o<^<i, ^ = ion%,4 ^r^l^' ^r^^2' 

where 3/4 < 5, p < 1. 

By c/>o we denote the cut-off function associated with the integral domain - the ball centred at x = 
of radius Rq - satisfying the above properties. 

Comparisons will be necessary between averaged quantities localized to cover elements at some 
scale R and averaged quantities at the integral scale. To accommodate this we impose several 
additional conditions for xi near the boundary of B(0,Rq). If B(xi,R) C -B(0, i?o) we assume 
ip < ipQ- Alternatively when B(xi, R) <f_ B(0, Rq) we stipulate that V = 1 on B(x , R) n B(0, Rq), 
satisfies (O, and we additionally have, 

ip = ipQ on the intersection of S(xq, Ro, 2Rq) and the cone with apex at the origin 
and with boundaries passing through the intersection of the circle centered at the 
origin of radius -Ro an d the boundary of B(xi, R), 

and, 

ip = on the intersection of the three sets B(0, R ) \ B{xi, 2R), S(0, R , 2R ), and 
the outside of the cone with apex at the origin and boundaries passing through the 
intersection of the circle centered at the origin of radius Rq and the boundary of 

B( Xl ,2R). 

The above apparatus is employed to study properties of a physical density at a physical scale R 
associated with the integral domain B(0,Rq) in a manner which we now illustrate. Let 6 be a 



4 



Z. BRADSHAW AND Z. GRUJlC 



physical density (e.g. a flux density) and define its localized spatio-temporal average on a cover 
element at scale R around 



1x t ,R=h [ 4 I 6{x,t)4(x,t) dx dt, 

1 JO -tl JB(xi,2R) 



lB(xi,2R) 

where < 5 < 1, and let (@}r denote the ensemble average over localized averages associated 
with cover elements, 

1 n 

(&)r = -J2°^r- 

8=1 

Examining the values obtained by ensemble averaging the averages associated to a variety of 
covers at a fixed scale allows us to draw conclusions about the flux density 6 at comparable and 
greater scales. For instance, stability (i.e. near constancy) of {(6) r} indicates that the sign of 6 is 
essentially uniform at scales comparable to or greater than R. On the other hand, if the sign were 
not essentially uniform at scale R, particular covers could be arranged to enhance negative and 
positive regions and thus give a wide range of sign varying values in {(©) r}. In order, then, to 
show the essential positivity of an a priori sign varying density 9 at a scale R it is sufficient to show 
the positivity of {(©)r}- 

Conversely, if 6 is an a priori non-negative density, then ensemble averages are all comparable to 
the integral scale average across the range < R < Rq. Precisely put, there exists depending 
only on K\ and K% so that, 

(3) ^6 < (Q) R < K-,00, 



where, 



1 f T 1 



9 = -/ ^3/ 6{x,t)0° o (x,t)dxdt. 



lB(0,2Ro) 

The inequalities (O follow in particular from the selection of our cut-off functions as well as the 
defining properties of our cover. Indeed, by directly comparing the integrands of the involved 
localized spatio-temporal averages to the integral domain spatio-temporal average, one sees that 
that, 

(4) ^e < (&) R < k 2 q . 

For additional discussion of (K±, K2) -covers and ensemble averages, including some computa- 
tional illustrations of the process, see MDaGr3L 



3. Enstrophy Cascades. 



Our mathematical setting is that of weak solutions to the magnetohydrodynamic equations (cf . 
IISeTell for the essential theory). Define V = {/ G L 2 (R 3 ) : V • / = in the sense of distributions} 
and let V be the closure of V under the norm of the Sobolev space (H 1 (R 3 )) 3 . By a solution to 3D 
MHD we mean a weak (distributional) solution to the following coupled system (3D MHD), 

ii t -AM + (u. V) u - {b ■ V)b + VP = 0, 

b t -Ab + (u- V)6 - (b ■ V)u = 0, 

V • u = V • b = 0, 

u(x, 0) = uq(x) G V, 

b(x,0) = b (x) £ V, 
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where the magnetic diffusivity and kinetic viscosity have been set to one and P is the total pres- 
sure. 

Taking the curl of the above equations yields the following evolution equations for the vorticity 
and current, denoted oj and j respectively, 



(5) d t u -Au = (u- V)w - (oj ■ V)u - (b ■ V)j + (j ■ V)&, 



3 



(6) d t j - Aj = (u • V)j - (j ■ V)u - (b ■ V)oj + (oj ■ V)b + 2 ^ Vfy x Vwj. 

z=i 

In our study we substitute for the inward kinetic and magnetic enstrophy fluxes through the 
boundary of a ball, B = B(xq, 2R), 

— j —\u}\ 2 (u-n)da = — j (u ■ X7)ui ■ ui dx, 
JdB 2 Jb 



[ \\j\ 2 i u - n ) d(J = ~ [ (u-V)j-jdx, 

JdB 1 JB 



the inward kinetic and magnetic enstrophy flux through a shell, S(xq, R, 2R), by incorporating a 
(nearly radial) cut-off function, </>. This cut-off function was defined in Section 2 and is chosen so 
that the gradient is directed inward. After multiplying (u ■ V) u and (u ■ V) j respectively by <j> u 
and <{> b, we have the following realization of the local kinetic and magnetic enstrophy fluxes at scale R 
around the point x , 

: = / \\u\ 2 (u ■ V<f>) dx = - J{u ■ V)w • (<f>u) dx, 
®i ,R ■= f \\i\ 2 ( u ■ V0) dx = - !{u- V)j ■ (4>j) dx. 



Formulas for the localized enstrophy fluxes are realized via the non-linear terms (u ■ V) oj and 
(u • V) j by multiplying and © respectively by 4> oj and <\> j and integrating. In this manner we 
see that the localized kinetic enstrophy flux is given by, 

i^(t) : = J J i| w | 2 (n-V(/>) dx ds = J ^\oj(x, t)\ 2 tp(x) dx + J J \V uj\ 2 (p dx ds 

t r i rt /• 

oj\ (d s <j) + A(f)) dx ds — I (uj ■ V)u • cjxjj dx ds 



>o J 2 jo j 

+ j J (b- V)j • (4>uj) dxds- J J (j ■ V)6 • {<puj) dx ds 

= j ^\cj(x,t)\ 2 tp(x) dx + / / \Vu\ 2 4>dx ds 



1 

2 



o 
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while the localized magnetic enstrophy flux is given by, 

^\j\ 2 {u-V4>)dxds= I i 



i / i|j| 2 («-V0)dxda= /^|j(x,t)| 2 V'(x) ( ix+ / / \Vj\ 2 4>dxds 



^\j\{ds4> + dx ds — j j (ui ■ V)6 • <pj dx ds 
I (b • V)j • (</>j) dxds- J J (j • V)tt • (#) 

3 



//(EV..XV6,). 



((/>j) ds 



(8) = J ^\j(x,t)\ 2 4>(x)dx + j \Vj\ 2 <j>dxds 

+ £T» + iVf + L J ' + N 3 2 + X. 

To establish the existence of the enstrophy cascades we will show that the ensemble averages of 
localized spatio-temporal averages of the above densities associated to an arbitrary (if i , i^-cover 
are positive and nearly constant across a range of scales. To achieve this several assumptions are 
in order. Before proceeding, we introduce several quantities. Fix a value so that, 

K, > max.{{K 1 K 2 ) ll2 ,ZK 2 /±,K 1 }, 

and set, 

a = AKpK 2 , 

where Kp is a constant which will be quantified in Section|5j 

Our assumptions are as follows. 

(Al) Hybrid Geometric/Smoothness Assumption. It is assumed for some threshold M > 
that we have, 

\u)(x + y,t) -u(x,t)\ < \u(x + y,t)\\y\2, 

2 

provided \y\ < 2Ro + R$, x in {|Vu| > M}, and u(x + y) ^ 0. 

Remark 2. This assumption is less satisfying than its analogue in the fluid case of [DaGr3|l 
where it was sufficient to assume the numerically and observationally motivated assump- 
tion of coherence of the direction of vorticity. While there is current to current and vortic- 
ity to vorticity alignment in 3D MHD, such a geometric consideration does not appear to 
sufficiently deplete the numerous nonlinearities (the fluid approach breaks down when 
examining the term Vui x Vfy), which explains our need for a different assumption. 

(A2) Modified Kraichnan-Type Scale. Let eo, Eq, and Pq denote the time averaged total energy, 
total enstrophy, and modified total palenstrophy at the integral scale. Precisely, 

5 









1 


eo = 


T , 


r 

Jo 


M 








1 


E = 


T , 


f 

Jo 










i 


Po = 


T , 


f 

Jo 





1 dx ds, 

2 2 



{\u\ 2 + \j\ 2 )dx ds, 



</>o(|Vw| 2 + |Vj| 2 ) dxds + ^J ±\u(x,T)\ 2 i/> (x) dx. 
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The modification of palenstrophy is due to the nature of the temporal cut-off; in addition, 
note that the cut-off's are modified for technical reasons. 



Set, 



«.= l£l, 



and, 



eo 
Po, 

Define the modified Kraichnan-type scale cto by, 

cr = max{4, e }. 

Our assumption (A2) is that, for some constant /3 = /3(M, ETi, K2, J Q T \\uj\\ 2 dt) with < f3 < 
1 (the precise value will be identified later), we have, 

Remark 3. The Kraichnan-type scale determines the lower limit of scales at which the 
enstrophy cascade occurs. For us this is realized by restricting to scales R with ctq/ f3 < R. 
In comparison to the analogous and identically named parameter in the 3D NSE case we 
here see a correction of f3 by a power of 1/2 necessitated by the emergence of energy-level 
quantities in Section HI 

(A3) Localization and Modulation - Because j Q T ds is an a priori bounded quantity (cf. 
BSeTel ), for a given constant Co > there exists Rq so that for any Ro < R£ we have 

T \ V 2 1 

2 • > - 1 



IMI 2 ds) < 

L 2 {B(0,2R +R%)) J k 

The localization assumption on Rq, the radius of the integral scale, is that, for Co = a, we 
have Rq < Rq. A precise (up to certain parameters) value for a will materialize in the 
proof. 

The modulation assumption imposes a restriction on the time evolution of the integral-scale 
kinetic and magnetic enstrophies across (0, T) consistent with our choice of the temporal 
cut-off. Precisely, 

\co(x, T)\ 2 ipo(x) dx > — sup / \oj(x, s)\ 2 ipo(x) dx, 
2 s J 

\j(x,T)\ 2 ip (x) dx > ^sup J \j(x,s)\ 2 ip (x) dx. 

Remark 4. Regarding localization, we have essentially introduced an upper bound on the 
range of scales across which the enstrophy cascade transpires; it is likely smaller than the 
scale at which the cascade initiates. This restriction is technical; a lemma to be included in 
a future paper will show that the near-constancy across a bounded range of scales extends 
to a range beyond that bound. A speculative statement of this lemma can be found in 
HDaGr3l . 

Using the terminologies of the above assumptions we are ready to state our main result which 
establishes the positivity and near-constancy of the kinetic and magnetic enstrophy fluxes across 
a range of physical scales. 
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Theorem 5. If a weak solution u, b of 3D MHD satisfies (A1)-(A3) on B(0, 2Rq) x (0, T) with initial data 
in V then, 

for all 2^ < R < Ro and K* dependent only on K\ and Ki- 

Remark 6. It will be plain that the localization estimates to be presented in the following section 
imply that (Al) alone guarantees smoothness over the spatio-temporal integral domain; hence, 
we are effectively concerned with the global-in-space (M 3 ) weak solutions that are smooth over 
the integral domain. However, since we do not impose any boundary conditions on the integral 
domain, the control over the 'smooth' norms is strictly local. 

Before continuing to the proof of Theorem [5] we observe that enstrophy flux locality - kinetic to 
kinetic, kinetic to magnetic, and magnetic to magnetic enstrophy exchanges are predominantly 
between scales of comparable size - is an immediate corollary. Discussion of this corollary and its 
precise statement is withheld until Section 6. 

To prove Theorem^ we will confine the localized (to a ball of radius R centered at xq) kinetic/ magnetic 
enstrophy flux between scale- and cover-independent multiples of the localized (to the integral 
domain) total-palenstrophy, P . In particular, for an appropriate fixed quantity we show, for 
any cover element centered at xq of scale R, that the space-time averaged inward fluxes of the 
magnetic and kinetic enstrophies are both nearly constant, ie, 

In order to achieve this certain bounds are needed. This is the subject of Section |4] The proof of 
Theorem|5]is presented in Section[5j 

4. Bounds. 

In this section each of the terms from ((7|) and are bounded by quantities which can be related to 
eo, Eq, and Po via the apparatus of ensemble averaging with respect to (Ki, K2) -covers at scale R. 
Throughout, we limit our consideration to a fixed ball of radius R and suppress the corresponding 
subscripts. We label various scalars by K, K e , Ke, and Kp and note these are dependent on K\, 
K2 and quantities determined by structural properties of 3D MHD. 

Bounds for the linear terms in 10 and © follow simply from properties of the spatial cut-off (see 

my. 

H:=H» + W < J ^(luf + UI 2 ) dxds. 

Before proceeding to bound the nonlinear terms a digression is necessary to introduce the kine- 
matic framework derived in ICo94| and adapted to MHD in IWuOOII and IHeXin|l . Recall that the 
deformation tensor of the velocity field can be decomposed as the sum of a symmetric component, 
the strain tensor of u, S, and a skew component, j x •. Put precisely, 

(9) V« = 5+^x. 

The operators in the above decomposition have the following singular integral representations: 

(10) u>(x) = ^ P.V. J a{y)uj{x + y)^, 
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and, 

(11) 

for, 



3 

4-7T 



P.V. J M(y,u(x + y)) 



dy 



-,a(y) = 3y®y- I, and M(y, f) = - (y <g> (y x /) + (y x /) <g> y) . 



A key feature for our treatment of the term Vui will follow from the fact that a and M (the latter 
when / is held constant as a function of y) have mean zero on the unit sphere. Integral operators 
such as these are the topic BSt93M and MSt70L We connect the above to the term Vu; x Vbi by noting 
for the unit vector e\ we have, 

Vu; = Vu e t + x e h 
Using the zero mean value property we write, 



|V^| < 



P.V. 

+ 

h + h 



1 3 \ dy 

— a(y) (u(x + y)~ U)(x)) + ^: M (y> ^(x + v) ~ ^( x ))j j^p 



y\<R 2 / 3 

[^-a(y)u(x + y) + ^-M(y,u;(x + y))) ^ 

\ y \>B?/ 3 V 47r 4?r / \y\ 



Our treatment is now divided between I\ and /2- For the former, the hybrid geometric-smoothness 
assumption (Al) entails that, 



h<K(a,M) [ \u( x + y)-u{x)\^ 

J\y\<R 2 / 3 \y\ 

dy 

'\y\<R?/3 



< K(a, M) i \u{x + y)[ 

J\v\<R 2 / 3 



y |5/2' 



and, therefore, by the Hardy-Little wood-Sobolev inequality (cf. Chapter V of |St93|), 



(12) 



/ills < K\\u\ 



L 2 (B(xi,R 2 / 3 ))- 



Regarding I2, Holder's inequality allows that, 

1 u{x + y) 



h < K(<j, M) [ 
JW 



\y\>R 2 / 3 \y\ \y\ 



dy 



(13) 



<K(a,M)( [ AdyY 

\J\y\>R 2 /3 |?/| 4 ) \J\y\>R 2 /3 W Z 

1 K 

^ K ( a ' M ) ^1/3+2/3 I M U 2 (R 3 ) = -jj\M\V 



u{x + y)\ 2 d \* 



Note that we can apply the exact same argument to u alone to obtain, 



(14) lj(x) = ]-P.V. [ 

4?r J\y\<R 2 /s \y\ 4vr 



dy 1 f dy 
a{y)(uj{x + y)- u(x))-—^ + — / a{y)u(x + y)-r-rz 

<R 2 /3 m 6 47T J\y\> R 2/3 \y\ A 



K , 



< ^ll^ll J L2( J B(x i , J R2/3) + -^-IMIl 2 
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We are now ready to establish bounds on the coupled non-linear terms. We begin with the most 
involved, the term involving Vuj x Vbi, that labelled X, as this will illustrate many of the compu- 
tational steps necessary for the other terms. We will show that, 

Q5) J Jcj)j- Vui x Vbidxds< ^Q su Pll^ill! + ^ j (t>\Vj\ 2 dxd^j 



We begin by splitting the spatial integral into the regions where | X7u\ > M and the complement. 
Considering the complement, 

r-T r r-T 



f [ <t>j- Vui x Vbi dxds< [ [ M\cj) 1/2 j\\0 1/2 Vbi\ dx ds 

Jo J\Vu\<M JO J\Vu\<M 



The second integral above can be expressed in terms of energy and palenstrophy level terms. 
Indeed, for k, h, and / distinct elements of {1, 2, 3}, 



[ [ ^^{dihf dxds = - [ [ tff-Hididibi dxds- [ [ d l <j) 2p ~ l b l d i b l dx ds 
Jo J Jo J Jo J 

= lo I (p2p ~ lb ^ dkjh ~ &hjk) dx ds + \ \l S did ^ 2p ~ lb2 i dx ds 

(16) - 2 [ I <f p ~ l \ b W SJ J\ dxds + \[ j W 

(17) < f (2a^<f) 2p - 3/2 \b\) (a-^cf) 1/2 \Vj\) dx ds 



^Wbl 2 dxds 



+ ihJ T I <t> Ap ~ 3 \ b \ 2 dxds - 



Applying Young's inequality to the first term of ((17)) gives us a final bound for the case when 

|Vu| < M, 

[ [ (pj -VuixVbidxds<- I [ <J)\Vj\ 2 dx ds + [ f 2p_1 |j| 2 dx ds 
Jo J\Vu\<M a Jo J R Jo J 

+^ r [^-Vdxds. 



R 4 Jo 

(18) < lQ S up||V;2j||2 + jT j cf>\Vj\ 2 dxds 
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Looking now at the regions of high kinetic intensity - i.e. | Vu| > M - we split into two cases using 
the decomposition for Vui, 



T f 

/ <pj ■ Vui x V6/ dx ds 

J\S/u\>M 



< [ T [ hl^VkW^jl dxds 

JO J\Vu\>M 

+ r i h^vbiu 1 ' 2 ^ dxd S . 

JO J\Vu\>M 

For the integral involving I\, the hybrid geometric/ smoothness assumption plus the localization 
and modulation assumptions will serve to minimize palenstrophy level terms. Applications of 
Holder's inequality and the bound ((12)) and subsequently the Gagliardo-Nirenberg inequality and 
Young's inequality yields an initial bound, 

fX r / rT \ h / rT 



f f hi^vkW^jidxdsKxf [ \M 2 LHB(x R2/3) dsY ( [ W^vkWlw^jWldsY 

Jo J\Vu\>M \Jo y y " ' J \Jo J 

< ^snp\\^j\\l + K £ llVi^Vb^ds^ 
(19) < -(IsnpW^jWl + K f T WV^^VhWlds 



a\2 

rT 







+ K J W^VVkWl ds 

where the constant a emerges from the localization assumption, (A3). We further decompose the 
last two terms. For the first, by expanding the integral, using the bound (Tl6l ), applying Young's 
inequality, and employing the properties of our cut-off functions, we have the following string of 
bounds: 

— [ T \\v4>^vb l \\ 2 2 ds = — [ T I (di^fidjbif dxds 

a Jo a Jo J 



rT 
aR 2 Jo 



^-^2 / / ^HdAfdxds 



< J-2 £ J ^i&iivji dxds+-^£ J ididi^im 2 dx ds 

<- r f <p\vj\ 2 dxd s +-^- r i cp 4p - 3 \b\ 2 dx ds. 

a Jo J art J J 



The last term of (Tl9l can be expressed in a similar fashion as the above but with several addi- 
tional steps. To begin, observe that by the product rule and integration by parts, we have the 
identity, 

||^VV6i||I= / Hdidjbi) 2 dx = J {d l d ] <p){d j b l ){d i b l ) dx + f <Xc^M(^W dx 

+ I (dj&idjkXdidik) dx+ ( {di(j))(djdjbi){dibi) dx. 
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We next bound each of the terms in the above identity. For the first, by Young's inequality we 
have, 

K 

a 



{didj(j))(djbi){dibi) dx < K J \didj<f)\\djbi\\dibi\ dx 
-wj ^'^MldM dx 



K 



ct> p - 1/2 \d j b l ^U 3 - 1 / 2 \d i b l \\ dx 



R 2 

Moving on, since b is divergence free and assuming k, h, and I are distinct, we have, 
K f K f 

— / (dj4>)(djbi)(didibi) dx = — {d j <p)(d j bi)(d h j k + d k j h ) dx 

K f 

< — / \dj(/)\\djbi\\dhjk\ dx 



< J (^^" 1/2 l^ 6 'l)(^l S ^'*l) dx 
" / ^ P_1 I V ^| 2 dx + \j ^l Vi ' 2 ^ 



Using Young's inequality and properties of the cut-off function, the remaining term is bounded 

as, 

K f K f K f 

— / 4>{djdjbi)(didibi) dx < — I <j)(didibi) 2 dx + — (f>(djdjbi) 2 dx 

< — [ <^|Vjf dx. 
a I 



Combining the above (and observing a > 1) gives a final bound for the term involving I\\ 

[ [ cf>j-VuixVbidxds<—(lswp\\^j\\l+f [ <f>\V j\ 2 dx ds] 

JO J\Vu\>M a \Z s J Q J J 

(20, ^fj ndxis 



Turning now to the non-singular case of our decomposition of |V«z|, we again use a to de- 
emphasize non-localization-apt palenstrophy level quantities. Its emergence is forced upon an 
application of Young's inequality with a reciprocal cost to an energy level quantity. Using the di- 
rect estimate l(T3)l on | J2 1 ar »d then applying Holder's inequality, Young's inequality, and the same 
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steps used to obtain and proceed from ((T7)l , we see that, 

/ 2 |/ /2 VMI/ /2 j| MldsjU ||^j|| 2 ||0iv6,|| 2 <i s y 



J\Vu\>M 



11 ,, , 1 ,,9 olK f T ,, , i , ., 
< --sup V 2 J 2 + / ^V6j 2^ 

a 2 a R Jo 

(21) < sup | + £ J <P\Vj\ 2 dxds 



Combining the bounds (TT8l ), d20l ), and ((2T|) establishes the bound (fl5l) and concludes our discussion 
of the term involving Vitj x V6j. 

Bounds for the remaining four critical order nonlinear terms, Nf, Nf, N(, and N%, are now at- 
tended. The processes for estimating Nf and N% are identical up to labeling. We only illustrate 
the latter. Splitting the space integral, when |Vu| < M, we have, 



2V| = [ [ (j ■ V)u ■ 4>j dxds<M ( 1103 j\ | dx ds. 
Jo J\Vu\<M Jo 

When | Vii| > M, noting by direct comparison that \ j\ < y/E\Vb\, we bound the quantity, 



T r- 

/ \Vu\\cf)^VEVb\\(f)^j\dxds, 

J\Vu\>M 

in an identical fashion to the unique-to-MHD term. The resultant final bounds are, 



N 2= I f (J ■ V)w • 4>j dxds< — (\ sup ||V)3 + / / <P\Vj\ 2 dx ds 





(22) + r / 02p_i|jf dx ds + ^ So i 04p_ * dx 

and, 

Nf = J J(oj- V)u ■ 4>u dx ds < — QsupUV'^Hl + y j <p\Vu]\ 2 dx ds 



(23) + § r [ ^-vp dxds+^ r [ <^- 3 ^ dx ds. 



The space-time integrals of N( and Nf also enjoy formally identical bounding procedures. The 
spatial integrals are still split into regions depending on | Vu\. When | Vu| < M we have the point- 
wise estimate |cj| < 5 1//2 M and we obtain a scaled version of prior bounds. When |Vit| > M we 
use the kinematic estimate, ((14)) , on \to\. The same familiar argument now shows that N( and Nf 
are bounded by the same dominating quantity appearing in (|22)l . 



Kp 

a 




+ 


K E [' 


R 2 Jo 


+ 


a 2 K e 


R A , 
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To summarize, 



2p_1 (M 2 + UI 2 ) ds 

Examining now L w and L J , we note, upon integrating by parts, a cancellation occurs when we 
consider their sum: 

L:=L w + L j = -[ f (j-tj)(V(j>-b)dxds. 



Although the above term is lower order and can be bounded in a more efficient way than what 
transpires below, we choose the less direct approach to limit problem specific dependencies of the 
parameter /3. We again split the spatial integral to obtain on one hand that, 

(j-uj)(V<j)-b)dxds<V5M- / ||026|| 2 ||^2j|| 2 dx ds 



J\Vu\<M R 



<^l !<? p - 1 \3\ 2 dxds + ^ I U^\b\ 2 dxds. 



R 2 Jo J 11 R A 
On the other hand, using the kinematic decomposition of to as a singular integral, we have, 



!! I a{y)u{x + y)^\j\\V<t>-b\dxds<-su V \\^j\\l 

Jo J\Vu\>M J\y\>R 2 / 3 \y\ a s 

+ j^- 3 \b\ 2 dxds, 



and, 



a{y)uj(x + y) -^r\j\\V(l)-b\dxds < -( / ||<^j|| 2 ||V(026)|| 2 dx ds 

J \ Vu\>M J \y\<R 2 / 3 \y\ a 



i 



_> 



< ^( s upll^ill! + ^ 11^ ViHl dx ds} 



a 2 K e f T ln _ 3 2 



o 



+ ^ I ir p -Wdxds. 



Combining the above we conclude that L is also bounded by the quantity given in (|22|) . 

To complete our discussion of bounds, observe that each of the enstrophy flux terms can also be 
bounded in a manner which will allow us to establish the independent cascades of both kinetic and 
magnetic ens trophies. These terms are lower order and their treatment is essentially the same as 
that for L. Indeed, the kinetic enstrophy flux, F w , ultimately satisfies the same bound as the final 
bound in d23t while the magnetic enstrophy flux, F 3 , satisfies that from (|22)) . 



5. Existence of the Enstrophy Cascades. 



In this section Theorem|5]is proven. We work in the context of an arbitrary (K±, K 2 )-cover {B(xi, i?)}™ = 
of -B(0, Ro) and assume the premises of Theorem [5] hold. 
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First we establish lower bounds for the following quantities which are independent of i (note that 
xo was arbitrary in the previous subsection and the bounds were independent of xq; our subscripts 
here indicate localization around X{ at scale R): 

rp rp 

J $% itR ds = j ]^\io{x,T)\ 2 ^{x)dx + J \Vu\ 2 fodxd8 + Hi + Ni + Li + Xi-F?, 

£ K,r ds = J \\H?,T)\ 2 il)i{x) dx + £ J \Vj\ 2 4>i dx ds + Hi + Ni + Li + Xi- Ff. 

We limit our discussion to the first of these as the arguments are identical. Observe that, 

Kp ( l i 

Hi + Ni + Li +Xi - Ff < SUp {Wt jW^iBix^R)) + Wi>iU\\tf(B(x i ,2R))) 



a \ 
+ 



* 

•T 



J j 0i(|Vi| 2 + |Va;| 2 ) dxds) 
+ ^ f T f ^\\u\ 2 + \j\ 2 )dxds 



The properties of (K±, K 2 ) -covers (see allow us, upon taking ensemble averages, to pass to a 
lower bound involving only integral scale quantities, 

1 f T 1 a \ >> 1 p ts Kp p ^E j, °?K2K e 

fj T^ Rds / R -l^ Po - K2 — Po -^ Eo -^^ eo - 

At this point we specify the value for a, 



and recall that, 
Consequently, 

and, therefore, 



a = AKpKl 



K, > maxj^!^) 172 , 3^2/4,^}. 



a ~ 4Ki ' 



1 f T J , \ . 3 „ K 2 ^ £ „ a 2 K 2 K, 

Recalling now (A2) we have, 



a 2 K 2 K e 4 2 



i? 4 



e < ^a z K 2 K e P Q , 



and, 



^E < P 2 K 2 K E P . 



Therefore, noting < (3 < 1, 
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The parameter modifying our Kraichnan-type scale, /3, is now chosen so that f3 2 is small enough 
to satisfy 

2p 2 K l K 2 (K E + a 2 K e )) < 1. 

Granted this, we obtain, 

Establishing the upper bound is, by properties of (K\, K 2 )-covers, immediate because we chose 
so that 4K* > 3K 2 . Indeed, 

(± £ ^'l-; .,, ds^ < K 2 P + 2P < 4K*P . 

Because the argument for (&)r is identical we conclude that, 



6. Flux Locality. 

As mentioned previously, we can immediately deduce locality of the flux from Theorem [5] Flux 
locality, in the context of turbulence phenomenology, refers to the fact that the transfer takes place 
primarily between comparable scales. This is realized in terms of flux by the proposition that the 
time averaged flux at scale R is well-correlated only with the time averaged fluxes at comparable 
scales. While flux locality is phenomenologically accepted in the hydrodynamic case, there has 
been some controversy about the locality in plasma turbulence (cf . MAIDisst AlEylO | for a locality 



result on the energy level, as well as a discussion and references on the topic). It is thus important 
to additionally examine the ratio of kinetic to magnetic enstrophy fluxes as a possible source of 
non-local enstrophy exchange. 

Define the time-averaged local magnetic and kinetic enstrophy fluxes associated to the cover ele- 
ment around the point follows, 



1 

TJo J 

i r T ri, 



K,r = £ / / >l 2 ( w • dx ds = ^K,* 



K,R = * I I Sbf(« • V0i) dx ds = R 3 & R . 



TJo J 2' 

Further, define the ensemble average over a (Ki, K 2 )-cover of these time-averaged fluxes as, 



1 n 



R, 

n 

i=i 
n 



(* i )R = zV*L R = R 3 (&)R- 



i=l 



Using the clear relationships between the spatio-temporal and ensemble averaged terms and the 
time and ensemble averaged terms one can use the bounds established in Theorem [5] to directly 
verify the following theorem (for which the proof is omitted). 
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Theorem 7. Let u and b satisfy the assumptions of Theorem\5\and let R and r be two scales in the range 
cr //3<r,R< R . Then, 



1 



■(- 



0^)r {^h (*")r 2 



Remark 8. Note that along the dyadic scale - r = 2 k R - the locality propagates exponentially. 
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